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Abstract 

We examine the singular behavior from the endpoint region x — > 1 in parton distributions unin- 
tegrated in both longitudinal and transverse momenta. We identify and regularize the singularities 
by using the subtraction method, and compare this with the cut-off regularization method. The 
counterterms for the distributions with subtractive regularization are given in coordinate space by 
compact all-order expressions in terms of eikonal-line operators. We carry out an explicit calcula- 
tion at one loop for the unintegrated quark distribution. We discuss the relation of the unintegrated 
parton distributions in subtractive regularization with the ordinary parton distributions. 
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Parton distributions unintegrated in both longitudinal and transverse momenta are used 
in QCD to analyze hadron scattering problems with multiple hard scales and to describe 
infrared- sensitive processes. See [l|, HI, [3| for reviews and references. These distributions 
represent less inclusive versions of the ordinary parton distributions. Accordingly, due to 
the lack of complete KLN [3] cancellation, they are affected by singularities from the region 
X — > 1 [l|, [HI, 0] corresponding to the exclusive phase-space boundary. 

In traditional applications these singularities are handled by placing a cut-off on the 
endpoint region. The cut-off can be implemented as the Minkowski-space angle obtained 
by moving away from the lightcone the eikonal line in the matrix element that defines the 
parton density, as in 0] and in the method of [6|, recently revisited in [sl, 0]. It can also be 
implemented in terms of the infrared cut-off associated with parton showering algorithms. 



11, 12 



as in Monte- Carlo event generators that make use of unintegrated densities [10 
But while the cut-off regularization is well-suited for leading order calculations, it makes it 
difficult to go beyond leading accuracy. Furthermore, with this method the connection with 
ordinary parton distributions and the lightcone limit are not so transparent. 

A more systematic approach is based on subtractive regularization. A formulation of 
the subtraction method, suitable for treating eikonal-line matrix elements to all orders, is 



given in [13j. In this approach the eikonal line attached to the field operator in the original 
matrix element remains lightlike, but the singularities are cancelled by counterterms provided 
by certain gauge-invariant eikonal correlators. The purpose of this paper is to study the 



unintegrated quark distribution using the method [13|. To this end we analyze the structure 
of the endpoint singularity in coordinate space. We carry out an explicit calculation at one 
loop. This allows us to identify the counterterms, and provides support for an all-order 
operator formula with subtractive regularization. We present the analysis for the quark 
distribution, as this contains the main aspects of the endpoint dynamics, but this treatment 
can also be given in the case of the gluon distribution. 

The analysis is given in terms of nonlocal operators. The techniques of [3] are used 
to make contact with the operator product expansion in local operators. Expressing the 
integral of unintegrated parton distributions in terms of ordinary distributions involves in 
general nontrivial coefficient functions, as discussed in [l5| for 0^ theory in dimension d = 6 



and in [16| for the x — > gluon density. The subtractive-method result that we find has the 
distinctive feature that the dependence on the regularization parameters introduced by the 
counterterms drops out in the integrated parton distribution. 

Let us first recall the basic behavior near the endpoint for fixed transverse momentum 
k± and lightcone momentum fraction x [l|. The unintegrated quark distribution f{x,k±), 
computed in a quark target with an infrared regulator p [H, 0], has the one- loop form 

f(i){x,k±) = PR{x,k^) - 5{1 - x)S{k^) j dx' dk'j^ Pn{x',k'^) , (1) 



where 



Pr{x, k±) as 



1 



+ {regular at x — 1} 



(2) 



1 — X k]_ + p'^ 

The X — > 1 singularity is the endpoint singularity, and is present for any k±. A physical 
observable O is constructed by integrating / against a test function (f (specifying the final 



2 



state, hard subprocess, etc.), which yields 



O = J dx dk± f{i){x,k±) Lp{x,k±) 

= J dki_ [(p{x, k±) - (p{l, 0_l)] Pr{x, k^) . (3) 

While in the inclusive case, with ip independent of fc^, the x — > 1 behavior in Eq. ([3]) simply 
corresponds to the familiar 1/(1 — x)+ distribution from real/virtual cancellation, in the 
general case uncancelled divergences are expected from the endpoint region. 

We now proceed as follows. We compute this singularity in coordinate space, we apply 
the subtractive regularization method [13|], and then going back to momentum space we see 
the generalization of Eq. ([3]) associated with it. 

We begin by considering the matrix element 

f^y) = \Y. (P,^l^(z/)K/W7+K)(n)^(0)b,s) . (4) 

s 

Here the ^/''s are the two quark fields evaluated at distance y = {0,y~,y±) with arbitrary 
y~ and y±, p and s are the momentum and spin of the target, taken to be a quark with 
pfi — (^p+^rn'^/{2p~^),0±), and V is the path-ordered exponential 



Vy{n) = V exp (^Qs j dr n^A^{y + t n)^ . 



(5) 



In Eq. (jS]) n is the direction of the eikonal line and A is the gauge field, = A^t", with t°- 
the color generators in the fundamental representation. The unintegrated quark distribution 
is obtained from the Fourier transform 

fix,k^) = J e-^^y-+^''--y- f{y) , (6) 

with d = 4 — 2e the space-time dimension. 

Let us expand the matrix element (jlj) to one loop. In Feynman gauge, the endpoint 
behavior results from graphs with gluons coupling to the eikonal line. Fig. [I](a) and (b). We 
calculate these contributions working in dimensional regularization with dimension d, and 
regulating the coUinear and soft divergences by keeping finite quark mass m and gluon mass 
A. We take n to be lightlike, n = (0, 1, Oj_), as is the case for ordinary parton distributions, 
and will give later the extension of the result to 7^ 0. Because we work in covariant gauge 
we need not consider the contribution from the gauge link at infinity. 

We expand the path-ordered exponential in Eq. (j4]) to first order in gA, and represent 
the gauge-field two-point correlator as 

rf« y ^ d,^{l) e^["(^^-A^+-)+^-(^.-01 . (7) 

Similarly, for quarks 



00 /• ^4 
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(a) 



(b) 



FIG. 1: One-loop contributions to the quark density from eikonal-line couplings. 



We switch to new integration variables f , a by setting 



a 



V a 



(9) 



with 0<t'<l,0<cr< oo. Then the integrals in the momenta i, q and r, a variables can 
be carried out explicitly for the graphs in Fig. [1] in terms of Bessel functions. We obtain 



/(a)+(fe) 



4d/2-2^d/2-l 
1 



V 



dv - 
1 - ^' 



^ip-yv 2(^/2-1 / P_ 
/i2 



X 



[(-2/2 + i£)^2Y/A- 



d/4-l 



d 



-\ \2-d/2 



(10) 



where K is the modified Bessel function, F is the Euler gamma function, is the dimensional- 
regularization scale, and we have defined 

p2 = (l_y)2^2^^A2 _ ^^^^ 

The singularity for w — 1 in the integrand of Eq. ffTOj) is the endpoint singularity. From the 
Fourier transform ([6]) we see that v has the meaning of plus momentum fraction (p+— 
where is the gluon's plus momentum. The y~ integral from Eq. ([6]) produces a 5{v — x) 
in the first term on the right hand side of Eq. ffTOj) and a 5(1 — x) in the second term, thus 
leading to the singularity structure for x — > 1 schematized in Eq. ([2]). Eq. ffTOj) shows in 
particular that the v 1 singularity is present even with finite A and m regulating the soft 
and collinear regions. 

We can see the relation of this result with ordinary parton distributions by expanding 
the answer (fTOj) i n powers of the distance y'^ from the lightcone. This step is analogous to 
the technique [ij] to analyze nonlocal string-like operators. We use the representation of K 



—y^ + ie 

oo 
k=l 



d/4-1 



2 2 \ 

y p^ 



-y2 + is)) 



,d 



2l-d/2 - ^) (p2)'i/2-2 



(12) 



2'^/2-3r(^ _ 2) (-^2^^5)2-^/2 



k=l 
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with 

. _ r(rf/2 - 1) _ r(3-d/2) 

A;! r(A; + d/2 - 1) ' k\ T{k + 3- d/2) " ^ ^ 



Inserting Eqs. (lT2j),(lT3l) in Eq. ([TO]) we get 



4^/2-2^^^/2-1^ Jo l-v + ie \^ J ' 2' > 

f> r(2 - d/2) V{d/2 - 1) y ^^/2+.-2._^2 2.. 

+ fc! 4fc r(A; + rf/2 - 1) ly/^J 

^ 4^^/2-2-^ nd/2 - 2) r(3 - d/2) y ... 2 2.2-./2-.4 .14. 

+ A;!r(A; + 3-d/2) ^ V ^^ ) J- U^J 

The expansion (fT4|) separates long-distance contributions in In(/i2/p2) and short-distance 
contributions in Ind/^/i^). The first hne of Eq. shows that the endpoint singularity 
f — > 1 cancels for ordinary parton distributions. At leading power the endpoint singularity 
is associated with the coefficient function in the second line of Eq. fll4l) . The higher order 
terms in the expansion are 0{y'^)^, with k > 1. 

Consider now ra^ 7^ in the matrix element (jl]). In this case the integrals in r and a are 
not elementary, and lead to formulas in terms of parabolic cylinder functions. The result 
can alternatively be given as the following integral representation, 

^ 7 o—i-Kd/4 rl roo poo 



X 



2 p ■ n V + — — ) p"*" + (1 — f ) m n 

JP-y p-*"^T^/(4cr) _ pip-yv -i(y-nT)^/{Aa) 



JO JO 

+ 



This can be used to study the two lightcone limits — *• and —>■ 0. In Eq. f|T5|) the 
behavior of the integrand at i; — * 1, resulting from the r integration over the eikonal line 
(see Eq. ([5])), is regularized by ^ 0. This is precisely the method jl, [sl to handle the 
endpoint singularity, giving a cut-off in x at fixed k± of order 1 — x > k^/ ^/^^ where 
Tj = {p ■ nf' /v?. The parton distribution obeys renormalization-group evolution equations 
in the cut-off parameter j^, 0, [13] and depends on 77, also after integration over transverse 
momenta. 

The subtractive method [13], which we now apply, works differently. This is reviewed 
in [l|. The matrix element (jlj) is still evaluated at n lightlike. It is multiplied however by 
vacuum expectation values of eikonal lines, which provide counterterms for the subtraction 
of the endpoint singularity. The counterterms contain in general both lightlike and non- 
lightlike eikonals. For this reason we introduce the vector 

= (m+,m",Ox) (16) 
and the path-ordered exponentials V^(m), Vy{u), where V is given in Eq. (jSj) and y is the 
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lightcone projection y = {0,y ,0±). We consider the matrix element (Fig. [2]) 

^(subtr)/ \ _ \ ^ 

^^^24- {Q\Vy{u)V^{nmn)V^{um / my{u)V^ {nmn)V^ ' 

where n = (0,1, 0_l). The numerator in Eq. f|T7j) coincides with Eq. (jlj). The denominator 
is the subtraction factor designed to cancel the endpoint singularity. Below we verify the 
cancellation explicitly at one loop. The subtraction factor is constructed using the tech- 
nique 13| and depends on both the lightlike eikonal in direction n and the non-lightlike 
eikonal in the auxiliary direction u. The unintegrated quark distribution with subtractive 
regularization is obtained as the Fourier transform of the matrix element f[T7|) . This 
distribution depends on the direction u. However, the dependence on u cancels in Eq. ([Tj 
for yj_ = 0. 












u 




y 



FIG. 2: Coordinate-space matrix element for the quark distribution with subtraction factors. 

We now go back to momentum space. We evaluate the Fourier transform of Eq. f[T7j) . 
expanding to one loop. We introduce the regularization parameter (, defined by the supple- 
mentary eikonal in direction u, 

(-) 

At one loop the result from the matrix element in the numerator is of the form in Eq. ([1]), 
while the vacuum expectation values in the denominator contribute subtraction terms. Ex- 
plicit calculation gives 



Wn{x, k^,C) + Sik^) j dk'^Wnix, k'^,0 



(19) 



where, restoring now also the contributions of finite (1 — x) and including non-eikonal graphs, 
we have 



Pr{x, kj_] 



asCp r (1 - x)[{k']_ + 771^(1 - xf - 2xm^] 



27r2 



[k^ + im?{l — x] 



212 



2a;/(l -x) 
[k'j_ + m'^{l — x) 



(20) 
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and 

Here we have set A = 0, ci = 4. Note the endpoint singularity (1 — x)^^ x asCF/{'n''^k\) 
for X — *• 1 in Pr, and the corresponding subtraction term in Wr. The specific form of the 
counterterms in the second fine of Eq. f|T9|) comes from the subtraction factors in Fig. [2l 
In particular, terms in 5(1 — x) cancel between the vacuum expectation values in Eq. fll7l) . 
This refiects the fact that in the coordinate-space results f[TUl) . f[T5]) only the terms in e^^'^^ 
depend on y±, while those in e*^'^ do not. 

Observe that the unintegrated distribution in Eq. (HM depends on the parameter ( of 
Eq. (fT8|) . but upon integration in k± the ^-dependence cancels between the two terms in the 
second line of Eq. (fT9|) . This implementation of the subtraction method is to be contrasted 
with the cut-off method [sl, [8|, where a residual dependence on the cut-off parameter is left 
in the integrated distribution. 

We may now consider the analogue of Eq. ([3]) for a physical observable O. Similarly to 
Eq. ([3]), suppose integrating a test function (f over the distribution /(™btr) Using Eq. ([T^. 
we obtain 

O = J dx dk± fli^^^'^\x, k±} ip{x,k±) 

= Jdx dk^ {{Pn - Wn) ^{x, k^) - Pr <^(1, 0^) + Wr ip{x, O^)} 

= j dx dk^ {Pr [ip{x, Ox) - <^(1, Ox)] + {Pr - Wr) [ip{x, k^) - ^{x, 0^)] } . (22) 

Unlike Eq. ([3]), the endpoint behavior in Eq. fl22l) is regularized. For x ^ 1 the first term 
in the last line of Eq. (122|) corresponds to the 1/(1 — x)_|_ distribution, while in the second 
term the endpoint singularity in Pr is cancelled by Wr. 

In conclusion, we verify at one loop that the subtractive method, implemented in Eq. fITTI) . 
provides a well-prescribed technique to identify and regularize the endpoint singularities 
in unintegrated parton distributions. This is accomplished while keeping the eikonal in 
direction n exactly lightlike, in contrast with approaches based on cut-off regularization 0, 
0, S [9I, and canceling the singularities instead by multiplicative, gauge-invariant factors. 
The one- loop counterterms in Eqs. (fT9|) . (122!) . generated from the expansion of these factors, 
correspond to an extension for /c^ 7^ of the plus-distribution regularization, characteristic of 
the inclusive case. In general, the operator expression in coordinate space given by Eq. (fT7|) 
can be used to any order. 

The subtractive method provides an alternative to the cut-off method, most commonly 
used in this context. In the cut-off method the eikonal n is moved away from the lightcone, 
as in Eq. f|T5|) above, and the cut-off is given hy t] = [p ■ nf' / rt^ . The lightcone limit is not 
smooth, so that one does not simply recover the ordinary parton distribution by integrating 
over k±_ and taking 77 — > 00 [6|. This behavior can be compared with the dependence on the 
"gauge-invariant cut-ofF' parameter ffTSjl in the case of the subtractive method. As noted 
below Eq. f|T7j) and below Eq. fl2Tl) . the dependence on the non-lightlike eikonal, introduced 
in Eq. f|T6l) to regularize the endpoint at unintegrated level, drops out of the distribution 
integrated over transverse momenta, where such regularization occurs independently by KLN 
cancellation. 
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Different methods of regularizing the endpoint singularities in the unintegrated parton 
distributions will result in different coefficient functions for the y'^ expansion of the kind in 
Eq. (fT4l) . The observed cancellation of the C dependence in the integrated density could 
be seen as corresponding to a particular scheme choice, which may be advantageous for 
applications such as the construction of event-generation methods 15, 18, 1^. It can also 
be helpful for studying the re-expansion of unintegrated distributions in terms of ordinary 



distributions |14| . Il5l . Il6[ |. The possibility of constructing one such scheme appears to be a 



distinctive feature of the subtractive method compared to the cut-off method. 
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